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超 深井 下 射 孔 暗 侵 彻 超 强 砂 着 的 ALE 仿真 


窦 益 华 , 徐 浩 , 李 明 飞 


(西安 石油 大 学 机 械 工 程 学 院 , 710065 西安 ) 


摘 要 : 随 着 完 钻井 深 不 断 增加 , 储 层 岩石 压 实 度 、 强 度 和 地 层 压 力也 不 断 增 加 , 需 采 用 超 强 射 孔 弹 
射 穿 污 染 带 , 超 强 弹 - 超 强 靶 -高 围 压 组 合 将 增加 侵 彻 分 析 的 难度 。 应 用 LS-DYNA 软件 ,结合 ALE 
算法 ,以 某 油 田间 下 射 孔 为 例 , 建 立 无 围 压 和 50 MPa 围 压 下 HS45-5 型 超 强 弹 - 枪 - 液 - 套 管 -砂岩 三 维 
模型 。 通 过 建立 非 反 人 射 边 界 条 件 , 消 除 爆 达 疲 对 射流 成 型 和 效果 的 干扰 ;通过 调整 关键 部 位 网 格 密 
度 、 优 化 不 同 接触 面 间 的 网 格 协调 ,减少 数值 畸变 与 网 格 数量 ,提高 分 析 精 度 , 降 低 机 时 。 研 究 表 
明 :13 hs 时 ,射流 能 量 峰值 达 60 kJ, 占 总 能 量 的 23% ,随后 急剧 衰减 ,850 hs 时 衰减 为 零 ;16 hs 时 ， 
射流 速度 峰值 达 6300 m/s;50 MPa 围 压 下 穿 深 为 560 mm , 较 无 围 压 下 降低 了 20% 。 通 过 本 研究 可 
为 并 下 超 强 弹 侵 彻 高 国 压 超 强 靶 的 仿真 分 析 提 供 参 考 。 
-关键 词 : 射 孔 侵 彻 ; 超 强 弹 ;高 围 压 ;LS-DYNA;ALE 算法 
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ALE simulation of ultra deep downhole perforating 
projectile penetrating super strong sandstone 


DOU Yihua, XU Hao ,LI Mingfei 


(School of Mechanical Engineering ,Xian Shiyou University ,710065 Xian ,China ) 


Abstract :With the increase of drilling depth ,the compaction ,strength and formation pressure of reservoir 
rock are also increasing,so it js necessary to use super-strong perforating projectile to penetrate the pollu- 
tion zone ,and the combination of super-strong projectile , super-strong target and high confining pressure 
will increase the difficulty of penetration analysis. In this paper, using LS-DYNA software , combined with 
ALE algorithm ,taking an oilfield downhole perforation as an example ,the HS45-5super-strong bullet-gun- 
liquid-casing-sandstone three-dimensional model is established without confining pressure and 50 MPa 
confining pressure. The interference of detonation wave on jet forming and effect is eliminated by establis- 
hing non-reflective boundary conditions. By adjusting the mesh density of key parts and optimizing the 
mesh coordination between different contact surfaces ,the numerical distortion can be reduced ,the analysis 
accuracy can be improved ,and the machine time can be reduced. The results show that the peak value of 


jet energy reaches 60 kJ at 13 ks ,accounting for 23% of the total energy ,and then decays sharply and rea- 
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ches zero at 850 js. At 16 Ms ,the peak velocity of jet reaches 6 300 m/s. The penetration depth under 50 


MPa confining pressure is 560 mm ,which is 20% lower than that without confining pressure. The research 


in this paper can provide a reference for the simulation analysis of underground super-strong projectile 


penetration into super-strong targets with high confining pressure. 


Key words :perforation penetration ; super projectile ;high confining pressure;LS-DYNA;ALE algorithm 


由 于 钻井 深度 的 增 大 ,深井 钻井 液 与 储 层 接触 
时 间 较 长 . 储 层 污染 较 重 ,从 而 导致 使 用 常规 里 孔 弹 
根本 无 法 达到 穿 次 要 求 , 因 此 需 采 用 大 装 药 量 的 射 


度 尚 不 清楚 ,可 能 会 导致 现 有 的 研究 结论 与 实际 井 
下 穿 深 存 在 较 大 差异 。 本 研究 针对 这 一 问题 进行 研 
究 ,建立 50 MPa 围 压 下 HS45-5 型 超 强 弹 - 枪 - 液 - 套 


孔 弹 进行 射 孔 。 除 此 之 外 ,目前 射 孔 方 面 研究 主要 
模拟 地 面 打靶 的 无 图 奈 工 况 , 超 深井 下 储 层 肝 压 党 
超过 50 MPa, 因 此 需 开展 超 强 弹 侵 彻 高 围 压 超 强 讽 
有 限 元 分 析 。 超 强 靶 ( 超 深 储 层 岩 石 ) 和 高 围 压 将 
增加 侵 彻 分 析 的 难度 。 

=s 现 有 文献 针对 超 强 弹 侵 彻 高 围 压 超 强 靶 分 析 较 
汰 > 文献 [1 ] 主要 研究 了 相同 装 药 量 下 不 同 装 药 结 
构 疝 射流 侵 彻 混凝土 靶 穿 深 的 影响 。 文 献 [2] 应 用 
LSEDYNA 软件 对 某 种 型 号 射 孔 弹 超 高 速 金属 的 身 
流 形成 以 及 作用 于 钢 靶 的 过 程 进 行 了 数值 模拟 。 文 
献 [3] 进 行 了 射 孔 弹 侵 彻 砂岩 靶 实 验 ,分 析 了 药 弄 
让 材料 、 壁 厚 、 开 口角 度 等 因素 对 射流 穿 深 的 影响 。 
区 扇 [4] 采 用 LS-DYNA 软件 ,固态 金属 药 型 日 转化 
成 汇 态 高 速 金属 流体 及 液态 高 速 金属 流体 与 混凝土 
介 颖 发 生 侵 彻 的 过 程 。 文 献 [5] 应 用 Workbench 软 
任 - 建 立 了 射 孔 段 管 柱 三 维 有 限 元 模型 ,分 析 了 管 柱 
的 动力 响应 规律 。 文 献 [6] 应 用 Workbench 建立 了 
某 现 格 射 孔 弹 侵 彻 圆柱 体 应 力 砂岩 的 有 限 元 模型 。 
分 怕 了 砂岩 应 力 分 别 为 0.50 MPa 时 射 孔 弹 侵 彻 砂 
崖 的 深度 。 文 献 [7] 采 用 LS-DYNA 软件 ,模拟 了 和 拢 
形 射 孔 弹射 流 的 形成 过 程 及 其 钢 彼 侵 彻 过 程 。 文 献 
[8-11] 用 LS-DYNA 软件 ,对 固态 金属 药 型 日 转化 成 
液态 高 速 金属 流体 及 液态 高 速 金属 流体 与 混凝土 介 
质 发 生 侵 彻 的 过 程 进行 了 有 限 元 模拟 。 文 献 [12] 
应 用 LS-DYNA 软件 对 聚 能 装 药 在 不 同 起 爆 方式 下 
的 射流 形成 过 程 进行 了 三 维 数值 模拟 。 文 献 [13-14] 
用 AUTODYN 对 油田 使 用 的 射 孔 弹 的 爆 又 压 震 药 型 
单 过 程 和 射流 侵 彻 套 管 过 程 进行 了 模拟 计算 。 文 
献 [15-17] 主要 研究 了 AIE 算法 在 相关 领域 的 应 
用 ,给 出 了 相关 的 研究 方法 和 技巧 ,为 本 研究 提供 借 
鉴 和 参考 。 

综 上 所 述 , 目 前 文献 对 于 超 深 井下 射 孔 爆 又 

侵 彻 深度 研究 很 不 完善 , 射 孔 侵 彻 研究 主要 集中 在 
井深 较 浅 的 工 况 下 ,实际 超 深井 下 射流 成 型 及 侵 深 


管 -砂岩 三 维 模型 ,应 用 LS-DYNA 软件 ,应 用 ALE 算 
法 进行 有 限 元 分 析 。 研 究 通过 以 下 3 种 方法 :中 对 
无 围 压 工 况 下 射 孔 弹 爆 艇 过 程 进行 了 模拟 ;通过 
在 后 处 理 中 定义 载荷 曲线 ,对 靶 板 施加 50 MPa 围 
压 , 模 拟 超 深 井下 工 况 ;@ 通 过 对 超 深井 下 射流 侵 彻 
过 程 的 数值 模拟 ,分 析 了 射 孔 爆笑 过程 中 能 量 转 换 、 
射流 头 部 速度 变化 。 通 过 本 研究 可 为 实际 井下 超 强 
弹 侵 彻 高 围 压 超 强 靶 的 仿真 分 析 提 供 方法 参考 。 


1 ALE 算法 描述 流 固 耦合 动力 学 
方程 


LS-DYNA 有 限 元 模拟 中 网 格 运动 是 通过 控制 
网 格 的 运动 速度 "或 控制 参考 材料 速度 w 来 实现 
的 , 即 


o(X) = Xi] (1) 
人 ux (XD) XK) az， 
? ot ot ot t[X] 
(2) 

Oox | _ 0% ox(X,t) oF 
WT 

ox(X,t) oxX I 

a x l= 7 (4) 


式 中 :外 ,t[x] 为 流体 坐标 x 固定 的 情况 下 对 流体 坐 
标的 时 间 导 数 ;x(X,) 为 某 时 刻 % 固定 情况 下 的 流 
体 坐 标 。 
结合 式 (1) \ 式 (2) 式 (4) ,得 到 ALE 描述 下 流 
固 耦 合 网 格 速度 公式 为 
v= -Co (5) 
应 用 LS-DYNA 软件 ,结合 本 节 ALE 算法 的 内 
核 , 即 可 实现 固态 金属 药 型 单 转化 成 液态 高 速 金 属 
流体 及 液态 高 速 金属 流体 与 混凝土 介质 发 生 侵 彻 的 
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2 有 限 元 模型 建立 及 动态 仿真 过 程 


2.1 模型 参数 及 材料 参数 确定 


为 了 研究 超 深井 下 HS45-5 射 孔 弹射 流 能 量 转 
换 、 射 流速 度 变 化 及 穿 深 ,建立 射 孔 弹 - 射 孔 枪 - 射 孔 


窄 着 华 ， 等 ; 超 深井 下 射 孔 弹 侵 彻 超 强 砂岩 的 ALE 仿 大 iNaX IVi 
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液 - 套 管 -砂岩 肢 模 型 。 模 型 参数 ;井深 5 500 m 、 射 孔 
段 深度 5 000 m, 地 层 压 力 梯 度 1.0, 射 孔 位 置 围 压 
50 MPa, HS45-5 射 孔 弹 , 炸 高 21 mm; 套 管 外 径 
178 mm, 壁 厚 12 mm; 射 孔 枪 外 径 102 mm, 壁 厚 
6 mm。 建 模 所 需 材 料 特性 参数 如 表 1 所 示 ， 包 括 套 
管 参数 、 射 孔 枪 参数 、 射 孔 液 参数 、 砂 岩 靶 参数 。 


表 1 仿真 模型 材料 特性 


Tab.1 Simulation model material properties 
模型 组 成 材料 状态 方程 强度 模型 失效 模型 密度 /(g .em 一) 

射 孔 枪 32CrMn4 None Kinematic 通用 准则 7.85 

套 管 TP140 Gruneisen Johnson-cook 通用 准则 7. 83 

射 孔 液 NULL Gruneisen NULL 通用 准则 1.0 

砂岩 SANDSTONE Polynomial Johoson-Holmquist Johnson-Holmquist 2.63 

空气 AIR Polynomial NULL 通用 准则 1.29 x10-3 

2.2 有 限 元 模型 建立 及 网 格 划分 到 射流 侵 彻 过 程 ,对 射流 成 型 和 穿 深 效果 有 所 干扰 。 
= 因此 在 建 模 过 程 中 对 模型 周围 空气 ,以 及 枪 管 、 射 孔 


2@31 有 限 元 网 格 划 分 

一 在 形成 聚 能 射流 和 侵 彻 射 孔 枪 和 套 管 的 过 程 中 
存在 大 变形 和 高 速 流动 , ALE 算法 ( Arbitrary La- 
gange-Euler) 可 以 很 好 地 解决 大 变形 和 流 固 耦 合 问 
题 5 因 此 空气 炸药、 药 型 黑 、 射 孔 液 均 采用 ALE 算 
法 ; 枪 管 、 套 管 . 弹 壳 、 砂 岩 采 用 Lagrange 算法 。 枪 管 
与 套 管 枪 管 与 砂岩 、 套 管 与 砂岩 之 间 的 接触 采用 侵 
衔接 触 类 型 。 根 据 模型 与 装 药 对 称 ,建立 1/4 模型 ， 
细 兹 射 孔 弹 及 其 侵 彻 轨迹 周围 的 网 格 ,模型 分 配 8 
个 单元 (包括 空气 单元 ) ,生成 102 862 个 网 格 ,如 图 
1 全 示 。 通 常情 况 下 ,建立 完整 射 孔 弹 模 型 计算 文 
件 较 大 ,计算 时 间 达 72 h, 且 容易 出 现 计 算 不 收敛 。 
本 研究 通过 建立 1/4 射 孔 弹 模型 ,调整 关键 部 位 网 
格 密度 ,优化 不 同 接触 面 间 的 网 格 协调 ,网 格 数量 可 
减少 50% ,计算 机 时 减少 40% 。 


图 1 射 孔 弹 - 枪 管 - 射 孔 液 - 套 管 -砂岩 有 限 元 模型 


Fig.1 Perforating bullet-barrel-perforating 


fluid-casing-sandstone finite element model 
2.2.2 有 限 元 边界 条 件 设 置 
由 于 炸药 爆炸 会 产生 爆笑 波 ,在 单 枚 射 孔 弹 数 
值 模拟 中 , 受 模型 尺寸 限制 ,爆笑 波 界面 反射 会 影响 


液 . 套 管 、 砂 岩 靶 施加 非 反射 边界 条 件 (non-reflecting 
poundary) , 消除 爆 麦 波 的 影响 ,如 图 2 所 示 。 


(b) 
管 . 射 孔 液 套 管 、 砂 岩 非 反射 边界 


Fig.2 Non-reflective boundaries of perforating projectile ,air， 


图 2 射 孔 弹 、 空 气 、 


barrel ,perforating fluid ,casing,and sandstone target 
2.2.3 初始 条 件 及 围 压 设置 

模型 建立 完成 之 后 , 先 在 ls-prepost 后 处 理 软件 
中 选择 需要 添加 围 压 的 面 ,选择 定义 负载 类 型 ,之 后 
在 uedit 软件 中 修改 输出 的 K 文件 ,添加 50 MPa 围 
压 ,关键 字 如 图 3, 其 中 DEFINE_CURVE 关键 字 为 
定义 一 条 负载 曲线 。 在 文件 编辑 器 uedit 中 分 别 定 
义 一 些 关 键 字 如 时 间 步 、. 计 算 时 间 、 沙 漏 能 控制 等 。 
其 中 时 间 步 设置 为 0.9 ps、 计 算 时 间 设 置 为 900 us， 
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开启 沙漏 能 避免 计算 出 现 负 体积 。 

*DEFINE CURVE 

1 2 1.000 1.000 0.000 0.000 
0.000000000000E+00 50.000000000000E-5 
1.0000000000000E+9 50.000000000000E-5 
*LOAD SEGMENT SET 
$ SSID LCID SF AT DT 

3 1 0 
xcontrol dynamic_relaxation 

250 0.100E-02 0.995 100 0.900 


图 3 50 MPa 围 压 关键 字 
Fig.3 Keyword of 50 MPa confining pressure 


*CONTROL TIMESTEP 
0000 0900000 0 000 0 0 1 
*CONTROL TERMINATION 
900.0 0 0.00000 0.00000 0.00000 

图 4 时 间 步 及 计算 时 间 设 置 


Fig.4 Time step and calculation time setting 


0 0.400E-01 1 


2,3_ 人 金属 射流 形成 及 侵 彻 过 程 


| 达 如 图 5 - 图 8 所 示 ,炸药 产生 的 巨大 冲击 波 在 册 

线 方 向 大 加 给 金属 药 型 罩 施 加 高 温 高 压 , 药 型 日 微 

轴 向 顶点 到 侧 向 底 边 依次 发 生变 形 ,在 轴线 方 

向 卫 发 生 激烈 磁 撞 ,在 尾 端 形 成 直径 较 大 的 杆 体 .前 

端 形 成 细 束 金属 射流 ,射流 依次 穿 过 射 孔 枪 , 侵 彻 套 

壬 草 筷 液 ,砂岩 彼 板 。 
GN 
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> 图 5 10 hs 射流 形成 
mm Fig.5 10 hs jet formation 
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习 6 16 hs 射流 侵 彻 枪 管 
Fig.6 16 hs jet penetrates the barrel 
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Time=24.998 


到 7 25 hs 射流 侵 彻 套 管 
Fig.7 25 hs jet penetrates the casing 
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Time=29 


图 8 29 ns 射流 侵 彻 砂岩 
Fig.8 The 29 hs jet penetrates the sandstone 


3 ” 射 孔 爆 壹 瞬间 及 侵 彻 过 程 能 量变 化 
规律 


爆 麦 产物 和 冲击 波 将 药 型 淖 向 中 间 压 垮 , 形 成 
高 速射 流 ,实现 爆 受 与 射流 能 量 转换 。 图 9 图 10 给 
出 了 射 孔 弹 爆笑 瞬间 炸药 能 量 、 药 型 第 能 量 随 时 间 
变化 曲线 。 如 图 9 所 示 ,炸药 能 量 最 大 值 约 为 60 kJ， 
炸药 能 量 于 15 ns 内 完全 释放 。 如 图 10 所 示 , 药 型 
音 能 量 于 10 ws 左右 增 至 最 大 ,最 大 值 为 14 kJ , 占 炸 
药 能 量 的 23% 。 由 此 可 知 , 有 23% 的 能 量 用 于 侵 彻 
枪 管 、 套 管 及 砂岩 。 


Qi 


0 20 40 60 80 100 
时 间 /us 


到 9 炸药 能 量变 化 曲线 


Fig.9 Explosive energy curve 


(10kDJ) 
© 


m 0.08 
缉 0.06 
到 0.04 
™ 0.02 
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Fig.10 Energy change curve of drug type mask 

图 11 为 前 50 ps 药 型 电能 量 随时 间 变 化 的 曲 
线 。13 hs 能 量 达 到 最 大 ,16 hs 时 曲线 出 现 第 1 个 
抛 点 ,此 时 射流 侵 彻 枪 管 并 消耗 能 量 ,但 能 量 损失 较 
少 ; 穿 过 枪 管 后 射流 侵 彻 枪 管 与 套 管 之 间 的 射 孔 液 ， 
射流 能 量 曲线 较为 平缓 ,基本 无 能 量 消耗 ;23 hs 时 
曲线 出 现 第 2 个 抛 点 ,此 时 射流 侵 彻 套 管 ,能量 较 穿 
过 枪 管 时 消耗 更 大 , 射 穿 套 管 后 进一步 扩 孔 持续 消 
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耗 能 量 ;27 ps 时 射流 开始 侵 彻 砂岩 ,可 以 看 出 此 时 
的 射流 能 量 开始 迅速 下 降 。 图 12 为 整个 侵 彻 过 程 
能 量 随时 间 变 化 曲线 ,从 图 12 中 可 以 看 出 ,射流 能 
量 在 刚 开 始 侵 彻 砂岩 的 210 ps 内 能 量 急 速 下 降 ; 
210 hs 之 后 射流 开始 出 现 断 裂 ,能 量 下 降 减 缓 ;于 
840 hs 左右 药 型 站 能 量 接近 于 零 ,表明 侵 彻 终止 。 
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图 12 射流 总 体能 量变 化 曲线 


Fig. 12 Overall jet energy change curve 
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>> 
45< 射 流速 度 变 化 规律 
5 


. 丑 射 流速 度 直接 影响 射流 最 终 穿 深 ,因此 需 对 侵 
彻 时 程 中 射流 速度 进行 分 析 。 图 13 为 射流 初 形成 
时 网 速度 云图 ,此 时 炸药 能 量 完全 释放 , 药 型 单 尖端 
被 压 震 ,形成 射流 。 由 图 13(a) 可 以 看 出 ,最 大 轴 向 
速度 出 现在 射流 尖端 ,最 小 出 现在 射流 的 尾部 ;由 图 
13(b) 可 以 看 出 最 大 径 向 速度 出 现在 射流 两 侧 。 图 
14 为 射流 断裂 时 射流 速度 云图 ,可 以 清楚 看 到 速度 
梯度 的 存在 ,射流 速度 最 大 集中 在 头 部 。 

图 15 .图 16 为 前 50 ns 射流 速度 .加 速度 曲线 。 
由 图 15 可 以 看 出 ,16 ps 时 ,射流 速度 最 大 达到 
6.3 km/s, 对 比 图 16 加 速度 曲线 ,射流 加 速度 为 零 ， 
16 hs 之 后 加 速度 均 不 大 于 零 。 因 此 ,16 ns 为 射流 
速度 最 大 时 刻 ,16 ps 之 后 射流 侵 彻 射 孔 枪 、 射 孔 液 、 
套 管 及 砂岩 靶 。 由 图 15 可 以 看 出 :16 ~ 18 us 时 , 身 
流 侵 彻 射 孔 枪 ,射流 速度 下 降 3% ,射流 加 速度 由 0 
变 为 负 值 ;18 ~ 22 hs, 射流 穿 过 射 孔 液 ,射流 速度 无 
明显 变化 ,射流 加 速度 为 0;22 ~ 25 hs, 射流 侵 彻 套 
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管 ,射流 速度 下 降 8% ,射流 加 速度 再 次 变 为 负 值 ， 
且 数值 大 于 16 ~18 ns 时 加 速度 数值 ;25 hs 之 后 射 
流 侵 彻 砂岩 靶 ,射流 速度 直线 减 小 ,加 速度 为 负 。 
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图 13 射流 速度 云图 
Fig. 13 Jet velocity cloud 
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图 14 射流 断裂 速度 云图 
Fig. 14 Jet breaking velocity cloud map 
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图 16 50 ns 射流 加 速度 变化 曲线 
Fig. 16 50 hs jet acceleration change curve 

图 17、 图 18 为 整个 侵 彻 过 程 射 流速 度 、 加 速度 
曲线 。 由 图 17 可 以 看 出 ,在 260 ns 左右 ,射流 速度 
下 降 斜 率 有 明显 变化 ,射流 速度 下 降 趋势 变 绥 。 图 
17 中 ,260 hs 时 射流 加 速度 数值 明显 变 小 ,初步 分 
析 为 此 时 射流 由 于 速度 梯度 存在 发 生 断 裂 ,只 有 射 
流 头 部 继续 侵 彻 ,其 余 射流 断 片 不 再 消耗 射流 能 量 。 
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840 hs 之 后 ,射流 速度 仍然 达到 300 m/s, 此 时 射流 
仍然 继续 侵 彻 靶 板 ,但 不 再 继续 向 前 侵 彻 ,而 是 属于 
扩 孔 阶段 ,直至 将 射流 能 量 消耗 至 尽 。 
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图 17 射流 速度 变化 曲线 


Fig.17 Jet velocity change curve 
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图 18 射流 加 速度 变化 曲线 


es Fig. 18 Jet acceleration change curve 


CN 图 19 为 无 围 压 和 50 MPa 围 压 工 况 下 选取 不 同 
时 间 点 位 移 所 绘制 的 位 移 时 间 曲线 。 从 图 19 中 可 
以 沽 出 :无 周 压 下 射流 穿 深 为 700 mm;50 MPa 围 压 
下 里 流 穿 深 为 560 mm, 对 比 无 围 压 下 穿 深 下 降 
20 吧 。 图 20 为 无 围 压 下 HS45-5 射 孔 弹 亲 侵 彻 砂岩 
试验 报告 单 。 由 图 可 知 , 无 转 压 下 HS45-4 射 孔 弹 穿 
深 为 690 mm, 穿 深 误 差 为 1. 4% ,在 计算 至 840 ms 
后 ,2 种 工 况 下 射流 穿 深 几 乎 不 再 有 变化 ,此 时 射流 
处 于 扩 孔 阶段 ,可 以 终止 计算 过 程 。 
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图 19 无 围 压 与 50 MPa 围 压 冉 流 穿 深 
Fig. 19 No confining pressure and 50 MPa confining 
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图 20 H45-$ 射 孔 试验 结果 报告 单 


Fig.20  H45-5 perforation test result report 


样品 编号 ; WSK094-2020(Q) 


| 穿孔 深度 到 


本 研究 应 用 LS-DYNA 软件 ,模拟 超 深 井下 工 
况 ,建立 1/4 完整 射 孔 侵 彻 的 三 维 有 限 元 模型 ,并 在 
砂岩 外 围 施加 50 MPa 于 压 , 应 用 ALE 算法 ,给 定 炸 
药 参数 ,模拟 爆 译 过 程 ,真正 实现 了 爆 胡 过 程 药 型 
四 固 流 转化 .射流 高 速 侵 彻 套 管 的 大 变形 和 流 固 磷 
合 仿 真 , 无 围 压 下 射流 穿 深 与 试验 穿 深 误 差 为 
1. 5% ,得 到 如 下 结论 。 

1) 通 过 设置 非 反 射 边 界 条 件 和 调整 关键 部 位 网 
格 密度 ,建立 完整 的 有 限 元 模型 ,消除 爆 钞 波 对 射流 
成 型 和 效果 的 干扰 ;减少 数值 畸变 ,提高 分 析 精 度 ， 
降低 机 时 。 

2)HS45-5 射 孔 弹 最 大 能 量 为 60 kJ, 药 型 单 能 量 
约 为 14 kJ, 炸药 能 量 转 化 率 为 23% ,射流 最 大 速度 
为 6300 m/s, 当 射流 速度 低 于 500 m/s, 射 流 处 于 扩 
孔 阶 段 ,不 再 继续 向 前 侵 彻 。 

3)50 MPa 围 压 下 模拟 穿 深 为 560 mm , 较 无 围 压 
下 穿 深 下 降 20% ,射流 在 600 js 之 后 侵 彻 速度 明显 
减 慢 ,在 840 ns 之 后 侵 彻 深度 不 再 变化 , 且 射 流 能 
量 趋 近 于 零 , 因 此 在 数值 模拟 中 终止 时 间 可 设置 为 
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900 us。 

4) 本 研究 通过 ALE 算法 模拟 射 孔 爆 禾 过 程 及 
侵 彻 深 度 , 计 算 结 果 与 试验 数据 误差 仅 为 1.5% , 计 
算 无 围 压 下 射 孔 弹 的 穿 深 合理 ,因此 可 以 用 ALE 算 
法 模拟 超 深 井下 (50 MPa 围 压 ) 射 孔 爆 受 以 及 射流 
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